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Abstract 

The three dimensional U(l) Lattice Gauge, in the weak coupling limit, 
is dual to a Discrete Gaussian model. We investigate this dual model and 
use it to calculate properties of the U(l) theory. We find that, because 
of the nature of the dual model, its advantages are outweighed by large 
autocorrelation times generated when the dual system becomes disordered. 



1 Introduction 



Duality transformations have provided a useful tool for investigating many the- 
ories both in the continuum and on the lattice. The term duality has been used 
to describe many types of transformation. A common property of these trans- 
formations is a mapping from the strong coupling region of one theory into the 
weak coupling region of its dual. Duality for lattice models was first developed 
by Kramers and Wannier for the two dimensional Ising model M and has since 
been extended to other more complex lattice models [0, 0]. Here we will use the 
term duality to refer to the extension of this Kramers and Wannier duality. 

This type of transformation is appealing from a practical point of view. Each 
model has its own set of advantages and disadvantages making simulations diffi- 
cult for some values of the coupling while they become relatively easy to simulate 
in other regions. 

By relating observables in the original theory to observables in its dual it 
becomes possible to use the dual model for simulations where the original theory 
encounters computational difficulties. This is particularly relevant in the weak 
coupling limit. Here the fields on the original lattice becomes ordered while the 
dual variables become increasingly disordered. 

In this paper we investigate a compact U(l) lattice gauge theory in three 
dimensions by performing simulations on its dual, the discrete Gaussian model. 
The duality transformation for this model was first performed by Kogut et.al. M 
and gives a model similar to the Solid on Solid models of surface transitions, but 
in three dimensions. 



We consider a compact U(l) gauge theory in 3 dimensions with dynamical vari- 
ables Ui }lt defined on links, where i labels a site on the lattice and \i specifies the 
direction. Each link M can be expressed in terms of an angle — 7r < 9i ifX < ir. 



2 U(l) and the Discrete Gaussian Model 
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The sum in equation (0) is over all sites % and directions 7 to generate all the 
plaquettes on the lattice. This gives a partition function: 

z = r nddj ne-pzav-^)) ( 5 ) 

Although this theory does not experience a phase transition, and thus has no 
continuum limit, it does possess an interesting topological structure. Through 
two duality like transformation, this partition function is transformed into the 
partition function for a coulomb gas of unbound magnetic monopoles plus a free 
photon field ||, |2], |3|, |!J. This identifies the instantons as magnetic monopoles 
which are exposed by writing the plaquette angle @ as 

6 a = Qn + 2im a (6) 

where ©□ is in the range — n < ©□ < n and n n = 0, ±1, ±2. This splits 9 n into 
a continuous part and a part which cannot be written as a gradient. One can 
therefore define a dual current 

jz-fi,p — ^^up n z,up (7) 

where z is a site on the dual lattice. The monopoles are then identified by writing 

^p,3z,p = m z (8) 
m z = 0,±1,±2,±3 (9) 

When m z 7^ there is a monopole/anti-monopole at the (dual) site and the 
current j Zjfl exposes the Dirac strings connecting monopole anti-monopole pairs. 
Dirac strings can also form closed loops without any monopoles. With periodic 
boundary conditions these Dirac strings can form closed loops which wrap around 
the boundary. The monopoles are physical objects which can be associated with 
cubes on the lattice however the Dirac strings are gauge dependent so their path 
can be altered through gauge transformations. Gauge transformations, however, 
cannot unwind loops which span the boundaries so these loops contribute to 
a winding number which is gauge invariant. The theory has only a confining 
'plasma' phase for all values of (3 although the monopoles are found to disappear 
from the system above j3 ~ 2.3 [0, [^]. 

2.1 The Duality Transformation 

The dual model, exposed by Kogut et.al. @, can be described by the partition 
function 

z = e -W£exp (£ln[4(/?)]) n^( e ^W A /3^,7p) (10) 
{k} \ ° / □ 

Where ki np is the dual variable, lk{0) is the modified Bessel function, N is the 
number of links on the lattice and we have dropped the overall factors of 2ti 
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from the delta functions. As we explicitly sum over all values of 9^ to generate 
equation ([K]) it implicitly includes all possible winding values. 

By associating the k's with links on the dual lattice, rather than plaquettes 
of the original lattice, equation ([T0| ) can be interpreted as a spin theory with the 
interactions on dual links given by 

k z ,a 2^ af 7P^*>7P (H) 

where z is the site dual to i. The constraints on these dual interactions then 
become 

e^aApkz^ = (12) 



2.2 Solving the constraints &; Boundary Conditions 

As equation (|l|) is just the curl of k, these constraints are automatically satisfied 
by writing k as the gradient of a scalar field. Thus we can express the fc Z)a 's in 
terms of a (dual) site variable <j) z 

k z , a = A a (j) z (13) 

where G Z. However we must be careful to include all the allowed k z , a configu- 
rations. While equation (|13|) guarantees the constraints are satisfied it does not 
take into account effects of the periodic boundary conditions. 

For an infinite lattice we can write the k za as a potential because of Stokes's 
Theorem. This works for any closed loops within the lattice. However if we have 
periodic boundary conditions a contour crossing the boundary an odd number of 
times forms a loop with no surface associated with it. Stokes's Theorem therefore 
doesn't apply. Writing k z ^ a as the gradient of forces this loop integral to be 
zero however the constraints flT2|) do not impose this restriction on k za . Therefore 



while equation ([13]) is sufficient within the lattice, at the boundary we need a more 
general prescription for fc Z)Q in order to generate all the allowed configurations. 
For the boundary links we can write k z ^ a as: 

k Zta = A a (p z + Q z , a {z) (14) 

Using the constraints (O) on Q gives 

e^A^A^ + Q z , a (z)) = (15) 

Thus 

KQzM = ( 16 ) 

for a 7^ v, so Q a cannot change across each plane perpendicular to a leaving 
only three extra degrees of freedom Q j , Qj and Qk ■ Now if we consider the line 
integral for a loop C which crosses the \x boundary once we get 



k^-dl = (17) 

c 

= Q, (18) 



3 



These extra parameters therefore allow the fields to differ by Q across the 
boundary without any cost to the action. This winding is the three dimensional 
extension of the 'step free energy' often introduced for two dimensional SOS 
models [|§. In two dimensions the cyclic boundary conditions are related to 
pinning the boundary of the surface. The extra 'step free energy', Q is introduced 
to remove this pinning constraint in the hope of improved statistics. 

The minimum requirement for generating all the interactions satisfying the 
constraints in (|T0|) with periodic boundary conditions is therefore 

= A/A + Q„ (19) 

for links on the boundary and 

K» = A /A (20) 

for all other points. 

As the step free energies Q are related to the boundary conditions they become 
less important as the lattice size is increased. When (3 is small the fields are 
predominantly at one level and, using the language of SOS models, the system 
is 'smooth'. Non-zero Q values therefore increase the energy for every link on 
the boundary and are subsequently suppressed by a factor proportional to the 
lattice size squared. At large (3 the system becomes 'rough' with a 'thickness' 
that increases with /3 and with the lattice size. In this region changing Q will, on 
average, decrease the energy for as many interactions as it increases it, provided Q 
is less than the thickness. Correlation functions are unaffected by these changes in 
Q and with only three winding parameters, compared with the L 3 site variables, 
fixing Q should not influence correlation function measurements, even for modest 
lattice sizes. We therefore set Q = below although it can be replaced at any 
time if desired. 

oo 

Z = e-W £ exp(£ln[/ (Ap , z) (/?)]) (21) 

{0} = -OO l d 

Where we have a sum over configurations of site variables <p z with an interaction 
defined on links k p = A p <j) z . 



2.3 The Weak Coupling limit 

By considering the weak coupling limit Kogut et.al were able to simplify (pi|) 
considerably. In this limit the partition function reduces to a discrete Gaussian 
model. 

HW))™-^ + P-HW*) (22) 

To improve this approximation we can systematically add higher order powers 
of n. This approximation is equivalent to using the periodic Gaussian or Villain 
action [|K| . The (3 on the right had side of equation ( p2"D can therefore be identified 
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with f3 v , the Villain beta. In his original paper Villain determined an approximate 
relationship between the cosine beta, (3 and (3 V . 



2 In 



(23) 



MP), 

Where Iq(P) and Ii(/3) are modified Bessel functions. This relationship is found 
through matching the first terms in the expansion of the cosine and Villain Boltz- 
mann factors and is therefore equivalent to adding extra terms in equation 
A good approximation to the partition function is therefore 



Z = $>xp £Mj 



/0 [|](A^) 2 | (24) 



2.4 Wilson loops 

In the weak coupling limit Wilson loops become correlation functions of the dual 
variables. Defining a loop W(J) where the current Jj iA1 defines the contour with 
J itfl given by: 

{+1 if link % — > i + /i is in C . 
— 1 if link i + /j — > i is in C. (25) 
otherwise 

The expectation value of this loop is then 

W(J) = — TT [ n e - /3 "(E l , D 1 - cos (^.°)+ i E liM ^.f J »^) (26) 

Z J —TT 

After the duality transformation, in the weak coupling limit, this become 

W {J) = -Y,e^ M ' S(J)) (27) 
Z M 

where in order to satisfy the new constraints is written as 

k^u = Sp,, + e^xAx4>i (28) 

Where Sn = 1 for an oriented surface with J p ^ as its boundary and zero 
everywhere else. Expanding the exponential for the extra terms due to Sn = 1 
leaves 

W(J) = We^(l+ +■■) (29) 

Z {<!>} V 2 Pv J 

where A is the number of plaquettes in the surface and the sum in each term 
is over the dual links l p = ^ep^S^ which penetrate the surface S^. Thus for 
the lxl Wilson loop, to order -4 

Pv 



W(D) = e~wz [\ + _ ^ _ (#m _ (30) 



where 4> = (</>). 
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3 Measurements 



In order to judge how useful the discrete Gaussian model is from a numerical 
perspective we also performed simulations using a U(l) lattice measuring the 
lxl Wilson loop directly. We considered lattices of size 8 3 and 16 3 with periodic 
boundary conditions. We used a standard metropolis algorithm and discarded 
the first 2,000 updates, where each update constituted a Metropolis hit for all 
links in the lattice. 50,000 updates were then performed taking measurements 
after every update. The simulations were performed on DEC Alpha 3400AXP 
workstations where, for a 16 3 lattice, each update took 0.212 ± 0.001 seconds. 

Before attempting accurate calculations for U(l) observables using the dis- 
crete Gaussian correlation functions we first checked that these measurements 
were insensitive to the step free energies. 

We compared simulations where the step free energy Q a was fixed at zero, with 
simulations where Q a was dynamically updated through a metropolis procedure 
before each measurements. We found that (|<5|) remained zero for /3 below ^2.2 
while above it increased monotonically. Measurements of the correlation functions 
however, were the same, within errorbars, for both static and dynamically varying 
step free energies. We therefore set Q a = for the remainder of our simulations. 

Simulations for the Discrete Gaussian model were performed on 8 3 and 16 3 
lattices. Again, in order to compare like with like, we implemented a metropolis 
algorithm with each update consisting of a metropolis hit on each site of the 
lattice. We ran 210,000 updates and discarded the first 10,000. Again we used 
DEC Alpha 3400AXP workstations and each update took 0.033 ± 0.001 seconds 
for a 16 3 lattice, almost 6.5 times faster than U(l). We measured the correlation 
functions after each update and used these to estimate the U(l) lxl Wilson 
loop through equation [3(]. The results are shown in figure [1] together with the 
strong and weak coupling expansions for the U(l) model []TTJ, |i~2]j . 

to order 0(— ) (31) 

P 

to order 0(/3 3 ) (32) 
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4 Errors and Efficiency 

In order to estimate the errors we split the raw data into bins of size 2 to 500. We 
used this blocking to measure the overall errors and calculate auto-correlation 
times r. We also checked that r calculated in this way was the same as T eX p 
calculated from an exponential fit to the auto-correlation functions. 

We used r to split the errors in the calculation into two parts, an autocor- 
relation (r) component and a standard error (cr) component. The total sample 
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Figure 1: (□) from U(l) and Discrete Gaussian simulations 



error, a s for a given simulation is then approximately 




where N is the sample size. The autocorrelation times measure the efficiency of 
the update algorithm and are shown in figure ^|a while a measures the statistical 
errors for each independent measurement and depends only on the model. These 
are shown in figure where again we have used blocking to estimated errors. 

For the U(l) lattice the measurement error decreases significantly above (3 ~ 
2. This is what we expect from the topological picture. At large f3 the monopoles 
condense out and the effective number of degrees of freedom therefore reduces 
considerably. Autocorrelations also increase around (3 ~ 2 signaling the emer- 
gence of extended objects. These are the monopoles which become important as 
the temperature decreases before they disappear from the system. Above (3 ~ 2 
the autocorrelation times fall when there are no monopoles. The autocorrelations 
remain short for large f3. 

For the Discrete Gaussian model autocorrelation times also increase around 
/3 ~ 2, however, they remain high as (3 increases. At small (3 the lattice is smooth 
but above f3 ~ 2 it becomes rough and the local algorithm slows down. This is 
different to the usual slowing down in spin models which occurs at low tempera- 
ture. Here the large scale objects of the dual model appear at hight temperature 
as the dual model becomes more disordered and the 'thickness' increases. This 
causes the local metropolis algorithm we used to have difficulty updating con- 
figurations. This suggest that a cluster algorithm may be effective at improving 
these simulations. A Solid On Solid algorithm has been developed for the two 
dimensional discrete Gaussian model by Evertz et al. This algorithm should 
be directly applicable to the three dimensional model although, to our knowledge, 
it has not been tested on a three dimensional lattice. The measurement error for 
the Discrete Gaussian model experiences a peak around (3 = 2 but most of this 
structure is generated by the transformation in equation (PU|). 

Figure [2] shows that both the auto-correlation times and the errors in each 
independent measurement are worse for the dual model for most values of beta. 
However the dual simulations ran over six times faster than the equivalent U(l) 
calculations so in order to compare the overall efficiency we need to consider 
simulations requiring the same amount of CPU time. Using equation we can 
scale the results to estimate the errors for simulations taking the same amount 
of computer time. This is shown in figure |3]. The dual model therefore only 
becomes more efficient for (3 < 2. However as the transformation utilizes the 
weak coupling limit the dual model is expected to be most accurate for large (3. 

5 Conclusion 

The main advantage of the dual model is the speed of its updates compared to the 
equivalent U(l) lattice. This increase in speed comes from using integer variables 
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Figure 3: Total error in (Sj/i) for a given CPU time. Errors are estimated using 
blocking 
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and a lack of local gauge invariance. This advantage however is destroyed by the 
large auto-correlation times generated by the local metropolis algorithm. The 
Discrete Gaussian model has the unusual characteristic of generating large-scale 
structures when the fields become disordered. This is the region where the U(l) 
system is at weak coupling and ordered and is when the transformation is most 
accurate. 

Cluster algorithms exist for both the U(l) and Discrete Gaussian models. And 
the Discrete Gaussian simulations would have benefited dramatically from an 
efficient algorithm. However our main concern was to study the transformation 
as a tool to improve the simulations. We were therefore careful to consider 
both theories on an equal footing. Here we have seen that by studying the dual 
model the characteristics of the simulations can be changes significantly. However 
whether this is of any practical use depends on the form of the dual model and 
the update algorithms available in each case. 

For U(l) it is unfortunate that the weak coupling region corresponds to the 
'rough' region of the dual model. This means that the disadvantages of both 
models occur when the duality transformation is most accurate. The opportu- 
nity to utilize better statistical properties of the dual theory therefore does not 
occur and the only gains produced are from the smaller configuration space. This 
is not always the case. Recently M.Zach et.al |14| have shown that in four dimen- 
sions, the U(l) dual model can be significantly more efficient than direct U(l) 
simulations. 



The recent progress in duality transformations for non-Abelian models flip 
allows for a similar study of more complex dual theories. These more realistic 
models transform to dual theories without their large scale excitations at high 
dual temperatures. They should therefore not encounter the simulation difficul- 
ties experienced by the discrete Gaussian model. For U(l) in three dimensions, 
simulations are relatively efficient and leave little opportunity for improvement. 
The increased complexity and difficulties experienced with non-abilian models 
allows more scope for improvement than this simple U(l) case. 
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